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ABSTRACT 

We estimate the conversion factor relating CO emission to H2 mass, acoj in five Local Group 
galaxies that span approximately an order of magnitude in metallicity — M 31, M 33, the Large 
Magellanic Cloud (LMC), NGC 6822, and the Small Magellanic Cloud (SMC). We model the dust 
mass along the line of sight from infrared (IR) emission and then solve for the aco that best allows 
a single gas-to-dust ratio ((5gdr) to describe each system. This approach remains sensitive to CO- 
dark envelopes H2 surrounding molecular clouds. In M 31, M 33, and the LMC we find aco ~ 3- 
9 M0 pc~^ (K km s~^)~^, consistent with the Milky Way value within the uncertainties. The two 
lowest metalhcity galaxies in our sample, NGC 6822 and the SMC (12 + log(0/H) « 8.2 and 8.0), 
exhibit a much higher aco- Our best estimates are 0^0*"^®*^^ ~ 30 M0 pc~^ (K km s^^)^^ and 
"co^ ~ 70 M© pc-2 (K km s^i)"^ These results are consistent with the conversion factor becoming 
a strong function of metallicity around 12 -I- log(0/H) ^ 8.4 — 8.2. We favor an interpretation where 
decreased dust-shielding leads to the dominance of CO-free envelopes around molecular clouds below 
this metallicity. 

Subject headings: Galaxies: ISM — (galaxies:) galaxies — (ISM:) dust, extinction — ISM: clouds — 
ISM: molecules 



1. INTRODUCTION 

In the local universe, stars form out of clouds made 
of m olecular hydrogen (H2) (e.g., iFukui fc Kawamural 
120101) . Understanding the processes that lead to star 
formation on large scales requires measuring the mass 
and distribution of this gas. Unfortunately, H2 lacks a 
dipole moment, most molecular gas is found under con- 
ditions too cold to excite quadrupole emission, and the 
high opacity of molecular clouds prevents UV absorption 
studies from probing the bulk of the gas. Estimates of a 
galaxy's H2 distribution therefore rely on indirect trac- 
ers, most commonly the lower rotational transitions of 
CO. The conversion between CO intensity and H2 abun- 
dance has been the topic of a great deal of investigation. 
In particular, the effect of metallicity and the local ra- 
diation field on the CO-to-H2 mass conversion factor, 
acoHj have been studied for more than two decades. 
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We work with «cOi the conversion from integrated CO in- 
tensity to mass of molecular gas. A linear scaling relates aco 



Infrared (IR) dust emission is a powerful tool to ad- 
dress this pro blem. Dust is observed to be well-mixed 
with gas fe.g.. lBohlin et al.lll978l: IBoulanger et al.lll996D 
and may be mapped by i ts emission at IR wavelengths 
(e.g.. ISchleeel et all [19981) . By compari ng IR emission, 
CO, and H I one can constrain aco fe.g.. lThr"on son 198^ 
[Israeli [r997b[ ). The procedure is to estimate the dust 
mass from IR emission, measure atomic gas (H I) and 
CO emission over a matched area, and then assume that 
the total gas mass (H2 -|- H I) is proportional to the dust 
mass, with the two related by a fixed gas-to-dust ratio. 
With measurements that span a range of relative CO, 
H I, and dust masses, it is possible to simultaneously con- 
strain the gas-to-dust ratio and aco- Key advantages of 
this approach are that it remains sensitive to any CO- 
free envelopes of H2 and that the calibration of aco is 
pinned to the H I within the system being studied. 

This exercise ha s bee n appli ed to several dwarf 
irregu l ar galaxies (llsrael Il997bl fal: iLerov et"aLl 120071 
120091: iGratier et al.l l20ld Bolatto et al. 2010, in 
preparation). In the low-metallicity Small Mag- 
ellanic Cloud, the results suggest a very large 
aco ~ 90-270 Mq pc'^ (K km s'^)-^. An analogous 
application to the Milky Wa y yields aco ~ 4- 
9 pc--^ (K km s-i)-i iBloemen et al.l [l990l: 

iDame et all[200l . roughly compatible with determina- 
tions of aco frorn fittin g the diffuse 7-ra y background 
([Strong fc Mattoxl[T996l : lAbdo et al.ll20Tol) . 

Most studies of aco have focused on a single galaxy or 
high l atitudes in the Milky Way, where confusion is min- 
imal. [Israeli (|1997bf) worked with a varied galaxy sam- 
ple but since his work the available IR and CO data 
for nearby galaxies have improved dramatically. This is 

to Xco, the conversion from integrated CO intensity to column 
density of H2. Including helium, Xqq cm^'^ (K km s~^)~^ = 
4.6 X lO^'Jaco Mo pc-2 (K km s"!)"!. 
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largely thanks to the Spitzer Space Telescope, which re- 
cently finished its cool mission and produced high quality 
maps of a number of Local Group galaxies. This is there- 
fore a natural time to apply this technique to measure 
aco across the Local Group in a self-consistent way. 

In this paper, we combine maps of CO, H I, and IR 
emission to estimate aco in five Local Group galaxies: 
the massive spiral M 31, the dwarf spiral M 33, the dwarf 
irregular NGC 6822, and the Large and Small Magel- 
lanic Clouds (the LMC and SMC). By treating aU five 
systems self-consistently, we minimize uncertainty in the 
" zero point" of the approach. With the Herschel mission 
now underway and the execution of several complemen- 
tary CO and H I surveys, this approach should be readily 
extensible to many nearby galaxies in the next few years. 

2. THE MODEL 

We and assume that dust and gas are linearly related 
by a gas-to-dust ratio, (5gdr, so that 



3GDR ^dust 



^H2 



^HI 



(1) 



where Sdust, 2]h2, and Shi are the mass surface densities 
of dust, H2, and H I along a line of sight. Substituting 
Sh2 = aco Ico we have 



^GDR Sdust = aco -^co + S 



HI 



(2) 



where the CO-to-H2 conversion factor, acOi and gas-to- 
dust ratio, (5gdr, are unknown and Sdust, Ico, and Shi 
are measured. After assembling Sdust, -^cOi and Ehi over 
many lines of sight in a region, we will use these data to 
solve for aco that best allows a single ^qdr to describe 
the data. 

2.1. aco-' Definitions and Scales 

The literature contains several working definitions of 
aco that apply to different scales or phases of the gas. 
In this paper and Equation [2] we define aco as the factor 
to convert from CO emission to total molecular gas mass 
on scales larger than individual clouds. Under this def- 
inition, aco includes any II2 associated with C^ in the 
outer, poorly shielded parts of clouds as well as gas imme- 
diately mixed with CO. Indeed, our goal is to measure 
whether such envelopes become dominant below some 
metallicity. Because this definition integrates over cloud 
structure, it is possible to derive a single aco for a whole 
galaxy or part of a galaxy, and to view that aco as a 
function of large-scale environmental factors. We choose 
this definition of aco because it is directly applicable 
to CO measurements of distant galaxies on kiloparsec 
scales. 

This is distinct from the ratio of CO to II2 only for 
molecular gas mixed with CO. Dynamical measurements 
using CO emission at high spatial resolution may mainly 
probe this quantity with limited sensitivity to H2 in an 
extended envelope not mixed with CO. The difference 
between this quantity and the quantity that we study 
has caused some confusion, leading to apparent contra- 
dictions between IR-based measurements and dynamical 
measurements. In fact these may be largely attributed 
to the different regions being probed (Section [B]). 

Similarly, we do not define aco as the ratio of H2 to 
CO along a pencil beam. We have no reason to expect 



TABLE 1 
Targets 



System Metallicity'' Resolution 
12-Hlog(0/H) (") (pc) 


M 31 








... inner 


9.0 


45 


170 


... north 


8.7 


45 


170 


... south 


8.7 


45 


170 


M 33 








...inner 


8.33 


45 


180 


...outer 


8.27 


45 


180 


LMC 


8.43 


240 


60 


NGC 6822 


8.2 


45 


110 


SMC 








... west 


8.02 


156 


45 


... east 


8.02 


156 


45 


... north 


8.02 


156 


45 



Ref erences: M 31 — |Y 
(I2OO9I) (compilation ); M 

IRosolowskv fc SimonI 

LMC and SMC— IDufouH 



et al.l 



33 

1 200a ): 
HID; 



Keller & Wo^ ((200^ : NGC 6822 
Israel (1997a) ( compilat i on). A dopted 
to match iBolatto et al.l ()2008l ) where 
possible. 

that this quantity, or any ratio that places many elements 
across a cloud, remains reasonably constant across part 
of a galaxy. To the contrary comparisons of dust and 
CO emission imply dramatic variations in the CO-to-Il2 
ratio within individual clouds (Pineda et al. 2008). 

Our spatial resolution ranges from 45-180 pc (Table 
[1]). For these resolutions, giant molecular clouds will 
mostly lie within one or two resolution elements (e.g., 
iHever et al.|[2009l) . This is ideal for our definition of aco • 
We wish to integrate over the structure of these clouds 
to make a single aco a more appropriate assumption. 

3. DATA 

To carry out this experiment we require maps of IR (to 
estimate Sdust), CO, and H I emission. Such maps have 
been pubhshed for M 31, M 33, NGC 6822, the LMC and 
the SMC. We refer to the original papers for details of 
the observations, reduction, and data. 

For M 31, we use the CO map taken bv iNieten et al.l 
(|2006[) using the IRAM 30 m telescope. This map 
has already been masked by INieten et al.l ()2006t ) and 
so contains only pos itive signal. We tra ce H I us- 
ing the 21cm map of Brinks fc Shand (jl984D . obtained 
with the West e rbork Synthesis Radio Telescope (WSRT) . 
iGordon et al.l ()2006D present Spitzer maps at 24, 70, and 
160/zm. 

For M 33, we take CO data from IRosolowskv et al.l 
(|2007[) . which combi nes BIMA (lEngargiola et al.l l2003D 
and FCRAO d ata (iHevCT et al.l 120041). We use the 
WSRT map bv IDeul fc van der Hulst (f987l). We use 



Spitzer data reduced following iGordon et al.l (l2005l) and 
presen ted by iVerlev et al.l ()2007[ ) and iTabatabaei et al.l 
(|2007l) . 

We use t he first CO map of the LMC obtained by 
NANTEN (jFukui et al.lll999( ). the Austraha Telescope 
Compact Ar ray (ATCA) -I- Parkes H I map of lKim et al.l 
(|1998l 120031). and IR maps from the Spitzer SAGE l egacy 
program " iMeixner et al.|[200l [Bernard et al.l [20081 ). 

We also use the NANTEN CO map of the SMC 
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(iMizuno et all 120011) . the ATCA + Parkes H I map by 
IStanimirovic et all (il999h . The SMC IR maps are a com- 



iGordon et alJ ( 


200f3, 


(IBolatto et alJ 


2007) 



and in prep.) and the S3MC survey 



For NGC 682 2 we use the IRAM 30m CO map 
bv IGratier et~all (|20H . the SINGS Spitzer maps pre- 
sented by Cannon et al.l (l2006f). and the VLA H I map of 
Ide Blok fc Walteii (|2003D . IGratier et"all (|2010f ) mapped 
the CO J = 2 — >■ 1 line for this galaxy, while th e rest of 
ourmaps are of CO J = 1 — >■ 0. We follow i Gratier et al.l 
(j2010( ) in assuming a line ratio of 0.7; we phrase all of 
our results in terms of CO J = 1 — )■ intensity assuming 
this ratio. We mask the CO map, keeping only emission 
above ~ 3cr at the original 15" resolution. 

For each galaxy we convolve all data to the resolution 
of the coarsest data set and align them on a common 
astrometric grid. In M 31, M 33, and NGC 6822 we are 
limited by the resolution of the Spitzer 160 /im data and 
so convolve all data to have a 45" (FWHM) Gaussian 
PSF. In the LMC the NANTEN CO data hmit our reso- 
lution to > 2'. 6; because signal-to-noise in the CO map is 
also a concern, we convolve all data to 4' resolution. The 
NANTEN CO data also set the resolution in the SMC, 
where we convolve all data to 2'. 6 resolution. 

In the LMC and SMC we subtra c t a fo reground from 
the IR maps. Following iBot et al.l ()2004D . we remove a 
scaled version of the Milky Way H I over these lin es of 
sight (for the exact approach see ILerov e t al. 2009). In 
M 33 and M 31 the reduction imposes the condition that 
the intensity is away from the galaxy, making a cirrus 
subtraction unnecessary. In NGC 6822, iCannon et all 
(|2006l ) already removed a Galactic foreground. 

The statistical uncertainty in the CO maps is about 
0.30 K km s-i in M 31, 0.35 K km s"! in M 33, 
0.30 K km s~i in the LMC, 0.01 K km s"^ in NGC 6822, 
and 0.08 K km in the SMC (though in each case 
this varies somewhat with position). The statistical noise 
in the H I maps is very roughly 1-3 Mq pc~^ with the 
uncertainty dominated by imperfect knowledge of the 
H I opacity and the reconstruction of extended emission. 
The noise in the IR maps is ~ 0.2 MJy sr~^ at 70 nm 
and ^ 0.7 MJy sr~^ at 160/um. The zero point in the 
IR maps is uncertain by ~ 0.1 MJy sr"-*^ at 70^m and 
- 0.5 MJy sr-i at 160/im. 

Throughout this paper IR intensity has units of MJy 
sr^^, color-corrected to the IRAS scale. HI surface 
density has units of Mq pc^^ and includes a factor of 
1.36 to account for helium. In the SMC a nd the LMC, 
H I in cludes an opacity correction based on lDickev et al.l 
(pOOOl ). In M 31, M 33, and NGC 6822 we have as- 
sumed that the H I is optically thin. CO J = 1 
intensity has units of K km s^^ and is related to 
II2 surface density, in units of M© pc~^, by a factor 
aco that includes helium (so that Eh2 = aco^co); 
if Xco = 2 X 10^0 cm-2 (K km s'^)''^ then 



aco 



4.4 Mq pc-^ (K km s^^)- 



4. METHOD 



Our goal is to identify regions where both H I and H2 
contribute significantly to the interstellar medium (ISM) , 
use the IR intensity in these regions to estimate the dust 
column, and then harness the assumption that gas and 
dust are linearly related to solve for aco- 



4.L Sampling and Target Region 

This experiment leverages our knowledge of H I column 
to infer aco- It thus works best in regions where both 
H I and II2 are important to the ISM mass budget. If we 
target areas where H I dominates then aco has little or 
no impact on ^gdr while if we target areas with only II2 
then aco and Squk are degenerate. 

For our targets and resolution, being dominated by 
II2 is not a concern for any reasonable aco- On the 
other hand each target has large areas where the ISM 
is overwhelmingly H I. Especially in M 31 and M 33, 
much of the H I is at large radius an d appears to have a 
different i^gdr from the inner galaxy (jNieten et al.ll2006D . 
We mostly avoid these H I-only regions, instead targeting 
the part of each galaxy near where CO is detected. This 
gives us a range of total gas surface densities and relative 
contributions by II2 and H I, allowing us to constrain 
aco and Jgdr- 

We define our target region by a ~ 3(7 intensity cut 
in the convolved CO map, Iqo > 1 K km s^^ in M 31, 
M 33, and the LMC, Ico > 0.25 K km s^^ in the SMC, 
and Ico > 0.03 K km s^^ in NGC 6822. We reject 
small regions (area < a resolution element) as likely noise 
spikes and then consider all area within about 1 resolu- 
tion element of the remaining emission. Finally, we re- 
quire a line of sight to have significant IR emission to 
be included; this is heo > 5 MJy sr"! in NGC 6822, 
I16O ^ 10 MJy sr"-'^ in all other targets. The result 
resembles loosely circling bright CO emission by hand. 
Figure [T] shows the target regions in contour on top of 
the 160/im map. 

In M 31, M 33, and the SMC treating the whole galaxy 
at once causes problems with our model, which assumes a 
constant aco and (5gdr across the ar ea studied. P r evious 
work arrived at similar conclusions. iNieten et al.l ()2006f ) 
observed (5gdr to be higher in the inner part of M 31 
than the 10 kpc ring containing most of the CO. In M 33, 
bright CO extends a fair distance out into the disk, but 
based on radial profiles of H I, CO, and 160/xm intensity 
it is immediately clear that this galaxy, too, has a strong 
radial gradient in 5gdr- The SMC's CO emission is clus- 
tered into three distinct regions that show evidence for 
local variations in their (^c nR.i aco. and giant molecular 
cloud (CMC) properties (|Lerov et al.ll2007tlMuller et al.l 
I20T0I) . 

To isolate regions with fixed aco and (5gdr, we sepa- 
rate each galaxy into several zones. We treat the "inner" 
part of M 31 separately from the 10 kpc ring and fur- 
ther divide the ring into a "north" and "south" part. 
We exclude 60° around the minor axis (in the plane of 
the galaxy) of M 31 (see We break M 33 into an 

"inner" zone where Tgai < 2 kpc and an "outer" zone 
where 2 kpc < Tgai < 4 kpc. We divide the SMC into 
three parts: a "west" region, a "north" region, and an 
"east" region. Our regions in the SMC deliberately ex- 
cl ude two clou d s nea r the center of the galaxy identified 
in ILerov et al.l (j2007l ) to have high CO-to-IR ratios; in- 
cluding these clouds leads to even worse solutions for 
this part of the SMC. Black lines and labels in Figure [T] 
indicate the divisions for each galaxy. 

After the target regions are defined, we sample each 
map using a hexagonal grid spaced by 1/2 the resolution 
(i.e., 22.5" in M 31, 20" in M 33, 2' in the LMC, 22.5" in 
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Fig. 1. — Images of our targets at 160/im (grayscale) with the region considered enclosed by a red contour. In M 31, M 33, and the SMC 
black lines separate regions that we treat separately and labels indicate how we refer to these in the text. 



NGC 6822, and 1'.3 in the SMC). This yields a matched, 
approximately Nyquist-sampled set of measurements of 
I70, I16O, Ico, and Shi- 

4.2. Estimating Dust Along the Line of Sight 

We follow the a pproach of IDraine fcTil ()2007D and 
IDraine et al.l (|2007| ) to estimate the amount of dust along 
the line of sight. These papers present models that can 
be used to estimate the dust mass and i ncident radi- 
ation field from IR intensities. Follow ing IDraine et al.l 
(|2007D and lMuiioz-Mateos etall ()2009D . we search a grid 
of models illuminated by different radiation fields to find 
the best- fit dust mass surface density, Sdust, for each set 
of IR intensities. A slight difference between our fits and 
those papers is that we quote the geometric mean of the 
dust mass across the region of parameter space where 
< Xmin + 1 within 1 of the for the best-fit 

model). We present the results of our own direct grid- 
search fits, but during anal ysis ma d e ex tensive use of the 
the work of Munoz-Mateo s et al.l (|2009), who parame- 
terized the results of model fitting over a wide range of 
parameter space have been as functions of the IR inten- 
sity at 24, 70, and 160^m. 

We only need a linear tracer of dust mass to solve for 
acOi the normalization affects (5gdr but not aco- While 
values of (5gdr that we find are interesting on thei r own, 
the overall normalization of the IDraine &: Lil (|2007| ) mod- 
els do not affect our resu lts for aco- Before settling 
on the IDraine &: Lil (|2007D models, we also used modi- 
fied blackbody fits with a range of emissivities to esti- 
mate the dust opacity, rigo- We obtained very similar 
results for acn using that approach, in the end prefer- 
ring the IDraine fc Lil (|2007l ) models mainly for their more 
straightforward handling of the 70/im band, which can 
include significant out-of-equilibrium emission. When we 
derive the uncertainties associated with our assumptions 
(Section 14.41 and Appendix |B]) we adopt either a mod- 



ified blackbody or the IDraine fc Lil ([2007*) models with 
equal probability. When using a modified blackbody fit, 
we include a variable fraction of 70/im emission from an 
out of equilibrium population and emissivity power law 
index in the calculation. 

M31: M 31 is quiescent compared to our other targets. 
It has a very low Ijq/Iiqo, implyi ng very low rad iation 
fields or colder dust te mperatures (|Gordon et al.l [20061. 
iMontalto et al.l (|2009l ) showed that with only Spitzer 
data the radiation field illuminating the dust is uncon- 
strained. In addition to poor constraints on the model, 
these low Iro/Iieo make ratio maps extremely sensitive 
to artifacts or contamination by point sources. M 31, of 
all our targets, will benefit most from the additional SED 
coverage offered by Herschel. Our approach in the mean- 
time is to treat all points in M 31 with as though they 
had the median IR color across th e whole galaxy. We 
implement the recommendation by IDraine et al.l ()2007f ) 
that in the absence of sub-millimeter data, the radiation 
field be limited to 0.7 times the local value. This treat- 
ment effectively assumes that value everywhere in M 31 
and uses the JigoMni to estimate the dust mass. We note 
this uncertainty in Table [5J 

NGC6822: NGC 6822 shows the faintest IR emission 
of any of our targets. As a result, many of the variations 
in /70//16O appear driv en by noise and arti facts rather 
than Tdust variations. iGratier et al.l (|201Cl( ) noted the 
difficulty of deriving dust temperatures for each line of 
sight. Based on plots of /70 vs. /leoi we identify two main 
" colors" in our data. We assign each line of sight the me- 
dian color for its group, with the cut at Ijo/heo — 0-6. 
We note this uncertainty in Tabled 

4.3. DijfuseRl 

Gas with different (5gdr may be superposed along the 
line of sight, in which case our assumption of a single 
<5gdr no longer applies. In some cases, we expect there 
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to be a diffuse H I component with httle or no associated 
dust along the line of sight. Such a component requires 
that we introduce an additional term into Equation [2] so 
that Sgas remains finite once Ejust 0. 

The simplest example of this is an edge-on spiral galaxy 
with a strong radial gradient in fc nR . M 31 is iiicline d 
and does show a gradient in 5gdr (jNieten et al.ll2Q06l ). 
The pile-up of many different radii along a single line 
of sight will be most severe along the minor axis and 
Equation [2] does not appear to describe these data as 
well as those along the major axis. Because we have 
plenty of data, we simply exclude 60° around the minor 
axes from our analysis. This improves the stability of 
our solution. No other correction appears necessary. 

The SMC has an extended distribution of H I that may 
also be very elongated along the line of sight. In this 
H I envelope and away from the main star-forming re- 
gions (il l the "Wing" and "Tail"), Jgdr is observed to 
be hi gh ()Bot et al.ll200l iLerov et al.ll2Q07l: iGordon et al.] 
I2009D . Similar results are found for othe r dwarf irregulars 
(jWalter et all 120071: iDraine et aLll2007[ ). Given the very 
high column densities found in the SMC it is plausible 
that low (5gdr H I lies along the line of sight. 

To account for an envelope of dust-poor H I, we sub- 
tract a diffuse component from the SMC H I map before 
fitting for aco ■ We estimate this component via an OLS 
bisector fit relating Sdust and Ehi at low Edust, where 
H I is likely to represent most of the gas. This line im- 
plies a value of Shi where Sdust — 0, which we take 
as our estimate of the diffuse H I. We estimate a dif- 
fuse component of Shi = 40, 20, and 60 Mq pc~^ for the 
western, northern, and eastern parts of the SMC. Similar 
fits to the other targets suggest a diffuse H I component 
with magnitude below 10 Mq pc~^, usually consistent 
with zero. Our uncertainty estimates include a 20% (Icr) 
uncertainty on this diffuse component in the SMC and 
±5 M0 pc~^ in the case of galaxies other th an the SMC. 

Thi s subtr action of diffuse H I differs from ILerov et al.l 
(|2007l 120091 though the latter applied this approach in 
some of the analysis). The effect on aco is largest in 
the "east" section of the SMC, which is embedded in the 
mostly diffuse Wing. 

4.4. Solution 

The Appendix lays out our method of solution and 
uncertainty estimates in detail. We identify the best- 
fit aco for each data set as the value that minimizes the 
point-to-point (RMS) scatter in log^o ^gdr- We estimate 
the uncertainty from three sources: 1) statistical noise, 
2) robustness to removal of individual data, and 3) as- 
sumptions. These are reported in Table [2] and we take 
the overall uncertainty to be the sum of all three terms 
in quadrature. 

4.5. Limitations of the Model 

Our model has several important limitations. It can- 
not recover a pervasive, CO-free II2 component like that 
suggested for the LMC bv .Bernard et al. (2008). To de- 
rive such a component we would need to make strong as- 
sumptions about (5gdr- We derive aco for H2 associated 
with CO emission on ~ 100 pc scales, essentially aco for 
GMCs and their envelopes. Based on UV spectroscopy 
(e.g., [Tumlinson et al...2002.) we consider a pervasive II2 



phase unlikely, but our experiment makes no test of this 
idea one way or the other. 

We assume that (5gdr does not vary between the 
atomic and molecular ISM. In fact, observations and the- 
ory suggest the (5gdr does correlate with density. The 
depletion of heavy eleme nts from the ga s phase increases 
with increasing density penkinsi 120091) , though absorp- 
tion measurements cannot probe to very high densities. 
Meanwhile, accounting for the observed dust abundance 
appears to require buildup of dust in molecular clouds 
fDwekl 11998: Zhukovska et al. 2008; Drainc 2009). A 
lower (5gdr in dense, molecular gas directly, linearly 
scales our derived aco- Dust already accounts for ~ 50% 
of the relevant heavy elements, so this effect cannot ex- 
ceed a factor of ~ 2 in magnitude. In fact, the ef- 
fect should be even less severe because we already study 
mainly the peaks of the gas distribution, where we expect 
H I to also have high densities. 

Our dust modeling also implicitly assumes that the 
dust emissivity does not vary between the atomic and 
molecular gas. Observations of suggest that dust prop- 
erties do change between the diffuse and dense ISM, with 
an enhanced emissivity in dense gas. The best evidence 
for this comes from a ^ 30-50% increase in far-IR opti- 
cal depth relative to optical extinction at high columns 
(^ 1 mag, e.g., |A rcc & GoodmanI 119991 : iDutra et al.l 
[200I iCamijresv et al.ii2005t) . Direct comparison of virial 
masses to s ubmillimeter em ission suggests a similar en- 
hancement (|Bot et al.ll2007| ) . The origin of this enhanced 
emissivity is thought to be a change in the size distri- 
bution of dust, with " fluffy," low- a l bedo grains created 
in dense environment s ()Dwek|[l997l : iStepnik et al.ll200"3l : 
iCambresv et al.|[2005() . As with Jgdr variations, emis- 
sivity variations directly scale the derived aco- 

Our best estimate is that emissivity and ^gdr varia- 
tions bias our measured aco high by a factor of ~ 1.5- 
2.0. Because this estimate distills a variety literature re- 
sults, none of them definitive, applying such a correction 
after the fact would simply confuse our results. More im- 
portantly, we have no handle on how these factors vary 
with environment, though we expect weaker effects at 
low metallicity, where shielding is weaker. We highlight 
a quantitative understanding of (5gdr and dust emissiv- 
ity vary with environment as the most important sys- 
tematics that must be addressed to improve dust-based 
derivations of aco - 

Finally, our data cannot rule out a pervasive popula- 
tion of very cold dust. We consider this unlikely and 
negle ct it throughout the paper (e.g., see IDraine et al.l 
I2007D . If such a population exists, the Herschel Space 
Telescope will identify it and our picture of dust will 
change dramatically over the next few years. However, 
we note that preliminary results suggest that dust masses 
are not dra matically affected b y the inclusion of Her- 
schel data ()Gordon et al.l [20T0I) . Other systematics in 
using dus t to t race II2 have been d iscussed at length 
elsewhere (jlsrae l 1997a. bl: [Schnee et"~ai . 2005; L erov et all 
120071 with references to many other works therein) and 
we do not repeat them here. Our best estimate is that 
none of these represent a major concern. Almost all will 
be addressed by the inclusion of long-wavelength Her- 
schel data (and are already improved relative to IRAS 
by using Spitzer^s 160/im band). 
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5. RESULTS 

Figure [5] plots scatter in logj^g ^gdr as a function of 
aco for our targets. In the background, a normalized 
histogram indicates the distribution of best-fit aco de- 
rived across the uncertainty exercises #2 (bootstrapping) 
and #3 (variation of assumptions). We find clear min- 
ima in each data set, which we report in Table [2] with 
associated uncertainties. We also quote the linear corre- 
lation coefficient between gas and dust and the scatter in 
logj^o ^GBR for the best- fit aco- Both quantities indicate 
the degree to which our model describes the system. 

The scatter plots in Figures [3l SI and [5] show gas {y- 
axis) as a function of dust (a;-axis) for each target. In the 
left column Egas is calculated from H I alone, the middle 
column shows Sgas = c^co-^co for our best-fit acOi and 
the right column shows total gas (Sgas = ^m + ctcolco)- 
In each case, the last column is a better match to a line 
through the origin (i.e., a single ^gdr) than the first two. 

The southern part of M 31 exemplifies the signal that 
we look for. H I correlates well with dust at low Sdust- 
At the high end, Sdust increases without a corresponding 
increase in H I (left panel). CO correlates well with Edust 
at high Sdust but drops to low values (or zero) while there 
is still significant dust, so that the intercept of the CO- 
Sdust relationship is not zero (middle panel). The gas 
becomes mostly molecular above a certain Sdust and is 
mostly atomic below this. Our best aco stitches these 
regimes together, allowing a single (5gdr to span all data. 

5.1. Comments on Solutions 

Before arriving at these results we varied our method- 
ology significantly. We tried different dust treatments, 
methods of solution, sampling regions, and treatments of 
diffuse H I. Most variations yield similar results to what 
we present here but a few cases are worth comment. 

First, the northern part of the SMC is not well- 
described by our model even after we excluded the high 
CO-to-IR clouds. Several distinct groups of data appear 
in the scatter plots for this region, suggesting local vari- 
ations in (^nn p or aon. The d ata for this region are 
already good (|Muller et al.ll2010l ). so improved modeling 
— multiple dust and Hi components, local variations in 
aco and Squr — seems like the most acute need. In 
the meantime our solution minimizes 5gdr and yields a 
result consistent with the other SMC regions. 

NGC 6822 yields reasonable solutions but is some- 
what unstable to choice of methodology. Adjusting the 
dust treatment, background subtraction, or weighting 
can change the best- fit aco by ~ 50%. This instabil- 
ity results mostly from the low signal to noise ratio and 
dynamic range in the data. These, in turn, are low be- 
cause of the large distance NGC 6822 compared to the 
SMC or the LMC, which makes observations challeng- 
ing. The resulting low intens ities also make con fusion 
with the Milky Way an issue () Cannon "eFall [20061) . Im- 
proved resolution from Herschel or ALMA should allow 
a greater range of Sdust as individual clouds are resolved. 
This will enable a stronger constraints on aco- 

The inner part of M 31 is also somewhat unstable to 
choice of methodology — especially region definition and 
fitting method. Best-fit solutions can range from aco ~ 
1-5. The most likely cause is variations in Jgdr and aco 
within the region studied — the data do not appear to be 



a plane in CO-H I-dust space. Both improved modeling 
and better IR SED coverage will help this case. 

We emphasize that our approach to the SMC is con- 
servative. We subtract a diffuse H I foreground, correct 
for high optical depth H I, and solve for (5gdr in the star- 
forming regions only — rendering us insensitive to any 
pervasive CO-free H2 component. All of these choices 
lower aco bringing it closer to the other galaxies, but 
the SMC still exhibits notably high a go. The mor e 
straightforward approaches employed in llsraell (|1997a[ ): 
iLeroy et al. (2007, 2009) yield even higher aco- 

5.2. aco ^'s. Metallicity 

Figure m shows aco (left) and Sgdr (right) as a func- 
tion of metallicity. Table [2] gives the adopted metal- 
licity for each target, which unfortunately represent a 
major source of uncertainty. Literature values span sev- 
eral tenths of a dex for most of our systems and internal 
variations in 12-|-log(0/H) often exist even where strong 
gradients are absent. We show a 0.2 dex uncertainty in 
12+log(0/H) in our plots to reflect this uncertainty. 

Figure [S] shows that (5gdr is a clear function of metal- 
licity, increasing with decreasing metallicity. This clear, 
one-to-one dependence of (5gdr on metallicity provides 
a consistency check on our solutions for aco- The 
data are roughly consistent with a linear relation (i.e., 
a flxed fraction of metals in dust), though a fit yields 
a slightly sublinear slope, logj^Q Jgdr = (9.4 ±1.1) — 
(0.85 ± 0.13) [12 + log(0/H)]. Our focus on star-forming 
regions likely biases us towards hi gh fcnR: studies of 
all material find a steeper relati on (jLisenfeld &: Ferraral 
ll998HMufioz-Mateos et alll2009f ). 

In the left panel of Figure |6] we plot the best- 
fit aco as a function metallicity. A gray area 
shows the c ommonly accepted range of values for the 
Milky Way (iBloemen et al.l 119861: iSolomon et all 119871: 
Strong fc Mattoxl 119961: iDame et all l2001t iHever et al.l 
20091). A dashe d line shows the aco argued by 
Draine et al.l (|2007[ ) to offer the best consistency between 
dust and gas measurements in SINGS galaxies. The 
molecular ring in M 31, both parts of M 33, and the 
LMC span a factor of ~ 2-3 in metallicity but show 
little clear gradient in aco- All three appear broadly 
consistent with the Milky Way. We find higher aco in 
NGC 6822 and the SMC, but the rise relative to the 
higher metallicity systems much steeper than linear. Fig- 
ure [5] therefore suggest a jump from "normal" to "high" 
aco at metallicity ^1/4 the solar value rather than a 
steady dependence of aco on metallicity. At the high 
metallicity end, the inner part of M 31 does suggest a 
low aco- It is unclear how much weight to ascribe to 
this point and whether the low aco is primarily a prod- 
uct of metallicity: the solution is somewhat unstable and 
a number of other conditions differ between the bulge of 
M 31 and our other targets. However, the point does 
highlight that lower-than-Galactic aco have often been 
observed, especially in ex treme, dust- rich systems (e.g., 
iDownes fc Solomonl[l998l ). 

6. DISCUSSION 

First, we emphasize a very basic result: combining H I, 
CO, and IR emission we obtain reasonable solutions for 
aco across the Local Group. In M 31, M 33, and the 
LMC aco appears consistent, within the uncertainties. 
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TABLE 2 

Dust-Based qco for Local Group Galaxies 



System 




Stat. 


Uncertainty'' 
Rob. Assump. 


Tot. 


Correlation'^ 


Scatter'' 
(logio ^GDr) 


Notes 




iVi ol 


















... inner 


2.1 


0.14 


0.03 0.26 


0.29 


0.68 ±0.02 


0.21 


U unconstrained 




... south 


10.0 


0.05 


0.01 0.12 


0.13 


0.85 ±0.02 


0.11 


U unconstrained 




... north 


4.2 


0.05 


0.02 0.19 


0.20 


0.83 ±0.02 


0.15 


U unconstrained 




M 33 


















...inner 


8.4 


0.07 


0.03 0.15 


0.17 


0.76 ±0.03 


0.08 






...outer 


5.0 


0.09 


0.05 0.17 


0.20 


0.76 ±0.03 


0.10 






LMC 


6.6 


0.04 


0.02 0.14 


0.14 


0.84 ±0.01 


0.14 






NGC 6822 


24 


0.13 


0.13 0.26 


0.31 


0.56 ±0.07 


0.18 


averaged colors 




SMC 


















... west 


67 


0.06 


0.02 0.08 


0.10 


0.91 ±0.04 


0.09 


subtracted Shi = 40 


Mq pc-2 e 


... east 


53 


0.08 


0.04 0.44 


0.45 


0.95 ±0.11 


0.07 


subtracted Sjji = 60 


M0 pc-2 « 


... north 


85 


0.16 


0.05 0.11 


0.20 


0.47 ±0.06 


0.15 


subtracted Shi = 20 


M0 pc-2 e 


Best fit. 


Units of M0 pc" 


-2 (K km s-i)-i 


including helium. 









^ Estimated la uncertainty in log^Q aco for each source of uncertainty. "Stat." — statistical (from Monte 
Carlo), "Rob." — robustness to removal of data (from bootstrapping), "Assump." — variation of assumptions. 

Linear correlation coefficient relating Egas and Edust for best-fit aco- Uncertainty gives la scatter for random 
re-pairings of data. 

la scatter in logjQ Sgdr in dex for best-fit aco- 
" Based on fit of Ehi to Edust at low Edust- 
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Fig. 2. — Scatter in log^o <5gdr as a function of aco; each dot shows a calculation for a trial aco- The shaded histograms in the 
background show the distribution of best-fit aco achieved across bootstrapping and variation of assumptions; they give some indication 
of the likelihood of a given aco describing the data. 



with the Milky Way. In most targets our aco agree 
well with previous determinations using other methods. 
Despite its simplicity and the potential importance of 
systematic effects, the approach outlined here produces 
consistent results. A lack of quantitative constraints on 
how the emissivity and Jgdr change between the atomic 
and molecular ISM (ideally measured at several metallic- 
ities) remains the most significant obstacle to derive ro- 
bust absolute values of aco from this approach. We cite 
circumstantial evidence that these effects combined may 
cause us to overestimate aco everywhere by a factor of 
~ 1.5-2 but emphasize the need for better constraints to 
improve the precision of this approach. In the meantime 
we have, we believe, an internally robust measurement 



of aco spanning an order of magnitude in metallicity. 



6.1. CO-dark Gas At Low Metallicities 

Does a rapid increase in aco over a narrow range 
in metallicity make physical sense? The metallicity 
dependence of aco is likely driven by "CO-dark" II2 
at low extinctions. In regions with low metallicity a 
significant mass of H2 is found in the outer parts of 
cloud s whe re the carbon is mostly C+ rather tha n CO 
(e.g. IMalon ev & Black 1988; Israel 1997b; Bolatto et al 
119991) Both .Glover k Mac Low, (,2010, ) and Wolfire et al 
(|201ClD recently studied this problem by modeling molec- 
ular clouds. They find that the dust extinction is the key 
parameter determining the fraction of CO-dark gas or 
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Fig. 3. — Scatter plots showing Sgas as a function of Sdust for M 31. The columns show the relationship between dust (ai-axis) and 
(left) H I, (middle) H2 derived from CO using our best-fit acOi a^nd ("S''*) total (H I + H2) gas. Our best-fit aco yields reasonable linear 
scalings between total gas and dust in each system (the red line shows the median (5qdj^), usually a clear improvement on the relationship 
between dust and either phase alone. 



aco- For clouds with mean extinction Ay ^ _8niag, such 
as those in the Milky Way, iWolfire et all (|2010f ) find a 
fraction of CO-dark gas foe ~ 30%, in ro ugh agreement 
with s everal obs ervations of nearby clouds (jGrenier et alj 
l2005t lAbdo et'a l. 2010). At intermediate metallicities 
the CO-dark H2 is not dominant, so that even a signifi- 
cant fractional change does not im pact agn very much . 
For example, in the calculations by IWolfire et al.l ()2010l ) 
going from the metallicity of M 31 to that of the LMC 
causes a doubling of the CO-dark gas fraction, from 
foG ~ 20% to foG ~ 40%, but the corresponding efi'ect 
on aco is only a factor of ~ 1.3, easily within the un- 
certainties of our study. As the mean extinction through 
the cloud decreases due to the effects of metallicity on 
the gas-to-dust ratio the increase in the fraction of CO- 
dark H2 makes it the dominant molecular component, at 
which point aco is very rapidly driven upward. 

Our measurements suggest that CO-dark gas becomes 



an important component in the metallicity range 12 + 
log(0/H) ~ 8.2 — 8.4. With such a smaU sample we 
cannot be sure that this is a general resul t, but it is in 
good agreement with the expectation from IWolfire et aTl 
(pOlQ) . Above these metallicities variations in the frac- 
tion of CO-dark gas will still exist, but their influence 
on the value of aco will be small. Excitation of the 
molecular gas (particularly due to temperature varia- 
tions) will be likely the dominant factor setting aco 
in molecule-rich systems. This is probably the cause 
of the low aco observed in LIRGs and ULIRGs (e.g., 
iDownes &: Solomonlll998») . which have solar o r slightly 
subsolar metallicities (e.g.. iRupke et al.] 120081 ). Within 
galaxies the radiation field incident on the molecular 
gas ma y also play a role, a fact emphasized by [Israeli 
(|1997al) . Unfortunately, the size scales involved make 
this ex tragalactic measurement difficult. While Israeli 
()1997a| ) found strong quantitative support for the infiu- 
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Fig. 4.— As Figure[3]for M 33 and the LMC. 



ence of the radiation fiel d on aco , iPineda et alJ ()2009l ) 
and iHughes et alJ ()201C1( ) recently failed to detect this 
effect in carefully controlled experiments in the LMC. 

6.2. Comparison With Literature 

Many authors have measured aco using a variety of 
techniques. We will not attempt to summarize the liter- 
ature, but focus on comparisons to two sets of measure- 
ments: 1) previous applications of IR-based techniques; 
and 2) high spatial resolution measurements of dynami- 
cal masses using CO emission. 

We measure high aco in the SMC and NGC 6822. 
This agrees with a larger trend in which IR photom- 
etry and spectroscopy suggest high aco in the range 
12-hlog(0/H) < 8.0-8.2. Israel ( 1997a) s aw this in a 
numb er of i rregulars. [ Rubio ct al. (200^, iLeroy et al.l 
(|2007D , and iLerov et"ai.i (,2009.) found similar results in 
the SMC, though as we note in Section [01 we actually 
solve for a somewhat lower aco than these studies — 
a fact we attribute to our focus on deriving Sqbr from 



the H I immediate associ ated with the star-forming re- 
gion. iGratier et al.l ()2010( 1 found the same high aco for 
NGC 6822 u sing several approach es. Using far-infrared 
spectroscopv iMadden et alT (|1997[ ) found indications of 
high aco in IC 10, a Local Group galaxy with metal- 
licity similar to NGC 6822 (we do not include IC 10 in 
this study because it lies near the Galacti c plane and has 
considerable IR foregrounds). iPak et all (|1998[ ) reached 
similar c onclusions in the SM C. D ust continuum mod- 
eling by [Bernard et al.l (|2008[ ) and iRoman-Duval et all 
(|2010D found a mixed picture in the LMC, also consis- 
tent with our findings. Note that in contrast with this 
broad agreemen t, recent work on diffuse lines of sight in 
the Milky Way (jLiszt et al.ll2010t ) suggests that the ratio 
of CO brightness to H2 column density is not a strong 
function of column density. 

Figure [7] compares aco as a function of metallicity 
between this study and the literature. The points in 
red indicate IR-based measurements. In detail, our mea- 
surements (circles) yield lower aco than previous IR- 
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based studies. We suspect that this is mainly because 
we solve for aco without assuming 5gdr or measuring 
it far away from the region of interest. One likely sense 
of systematic variations in Jgdr is that ^gdr is likely to 
be higher in the dense gas close to molecular complexes, 
which tend to reside mainly in the stellar disk, tha.n in a 
diffuse, extended H I disk fe.g.. iStanimirovic et al.lll999l : 
iDraine et afl 120071: iMufioz-Mateos et al.ll2009[) . If Jgdr 
is taken to be too high. Equation [2] yields a correspond- 
ing overestimate of acQ. In the SMC our attempt to 
remove a diffuse Hi component along the line of sight 
also leads to lower aco (Section 15. ip . though it is less 
clear that our approach is correct in that case. Regard- 
less of the cause, by simultaneously solving for aco and 
^GDR in the regions of interest in a uniform way across 
a heterogeneous sample we improve on literature studies 
of individual galaxies. 

A long standing discrepancy exists between IR-based 
results and high-resolution virial mass measu rements 
based o n CO observations . Using virial n iasses. I Wilson 
1995() . iRosolowskv et al.l (120031 ). and iBolatto et al 



20081 ) all found weak or absent trends in Xqo as a func- 
tion of metallicity. The blue points in Figure [7] show 
virial mass results from CO observations with resolu- 
tion better than 30 pc. The two approaches agree up 
to about the metallicity of M 33 or the LMC, and then 
strongly diverge in the SMC. This divergence is most 
easily understood if the additional H2 traced by IR lies 
in an extended envelope outside the main CO emitting 
region ([Bolatto et al.ii2008r) . Such an envelope can rec- 
oncile the virial mass and dust measurements and natu- 
rally explains the scale-dependence of aco observed by 
iRubio et all (|1993D in the SMC. These envelopes could 
perhaps still have an effect on the velocity dispersion of 
the material inside it (and consequently the measured 
virial mass) via surface pressure. Structures with virial 
parameters a < 1, however, are often observed inside lo- 
cal molecular clouds, suggesting that at least in some 
instances the velocity dispersion does not appreciably 
show the impact of the su rrounding material. A n al- 
ternative view is argued bv lBot et al.l ()2007l |2010), who 
observed discrepancies between dust-based masses and 
virial masses even at fairly small scales. They suggest 
that magnetic support becomes very strong at low metal- 
licities, perhaps because of higher ionization fractions in- 
side clouds. 



7. SUMMARY 

We combine CO, HI, and IR measurements to solve 
for the CO-to-H2 conversion factor, aco, in M 33, M 31, 
NGC 6822, the LMC and the SMC. We estimate the dust 
mass from IR intensities and then identify the aco that 
produces the best linear relation between total (H 
H2) gas and dust. We accomplish this finding the aco 
and (5gdr that minimize the scatter about Equation [5J 
We find that aco is approximately constant (within a 
factor of two) in M 31, M 33, and the LMC, with a 
value aco ~ 6 M0 pc~^ (K km s""'^)"^. By contrast 
NGC 6822 and the SMC, the lowest metahicity galax- 
ies in the sample, show a drastically higher aco, ~ 30 
and 70. The resulting gas-to-dust ratio, (5gdr, scales ap- 
proximately linear with metallicity. 

We attribute the behavior of aco to the transition 
from the regime where most H2 is bright in CO to a 
regime where CO is mostly photodissociated and the 
bulk of the molecular reservoir is CO-dark. In our 
sample, this transition occurs around 12 -1- log(0/H) ^ 
8.4 — 8.2. With only a limited number of systems, an ac- 
tual numerical prescription for aco is beyond the scope 
of the paper. These results agree qualitatively with a 
large body of existing work using IR based techniques, 
though quantitatively we find lower aco than previous 
work (e.g., Israel 1 997a ), probably because we restrict 
our analysis to CO-emitting regions. 
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APPENDIX 
METHOD OF SOLUTION 

For a given acoj we can compute the total gas surface density (Sgas — + Shi) for each hne of sight. With the 
corresponding Sdust, each point imphes a value of (5gdr- We have assumed (Equation [2]) that a single ^gdr describes 
each data set. A simple way to identify the aco that best fits this model is to try many aco and select the one 
that minimizes the scatter in logj^g (5gdr- First we pick a set of aco that bracket the true value. For each aco, we 
calculate Egas and combine it with Edust to measure the RMS scatter in logj^g (5gdr- This yields clear minima in each 
data set, which we identify as our best- fit aco- 

We experimented with several other approaches, including plane fits with errors in all three axes. The results agree 
with what we report here except that the plane fit is less stable in the inner part of M 31, NGC 6822, and the northern 
part of the SMC. The presence of significant intrinsic scatter in the relation makes the weighting of data subjective 
and so the plane fit is not significantly more rigorous than using the scatter as a goodness-of-fit metric. We also 
experimented with maximization of the linear and rank correlation coefficients relating dust and gas. Most of these 
approaches yield similar results. Our uncertainty estimates include substituting the median-absolute-deviation for the 
RMS in the goodness of fit. 

This method can fail in a data set without a good mix of H2 and H I-dominated lines of sight. For example, in a 
data set with 100 overwhelmingly H I lines of sight and 3 mostly molecular lines of sight, the H I-dominated data will 
drive the scatter in logj^Q i5gdr leaving little sensitivity to aco- In such a case we would like to weight the high H2 
lines of sight more heavily. In practice, our definition of sampling regions addresses this concern and we weight all 
points equally. 

UNCERTAINTIES 

We gauge the uncertainty in the best-fit aco in three ways: 1) We vary details of the calculation (the calibration 
and background for each map, the dust model, and the goodness of fit statistic) across their plausible range, repeating 
our solution for each new set of assumptions. 2) We solve for aco in ^ simulated noisy data set where we know the 
true aco by construction. 3) We repeat the original solution while bootstrapping, solving for aco using data drawn 
from the original sample (with the same number of elements) allowing repeats. We thus assess the sensitivity of aco 
to our assumptions, statistical uncertainty, and robustness. We bookkeep each as a gain-style uncertainty (i.e., RMS 
in log aco) and then take our overall uncertainty to be the sum of all three te rms in quadrature. 

In the first test, we carry out 100 it erations. Although the lDraine &: Lil ()2007[) models have been applied successfully 
to the SMC (jSandst rom et aTl I2010I ). adopting them represents a key assumption. Therefore, in each iteration we 
randomly select with equal probability to use either the Dra ine & L i (2007) dust models or a modified blackbody with 
/3 = 1-2 and assuming 50 % of the 70/xm ern ission to come from out-of-equilibrium very small grain emission, a value 
appropriate for the SMC (|Lerov et al.ll2009t ). We also select our goodness of fit statistic to be with equal probability 
either the RMS scatter (as in the main result) in log^g^GDR, the fractional scatter in (5gdr, or the median absolute 
deviation in logj^g ^gdr- We randomize the the zero point of the IR maps and the absolute calibration of each data set 
within their plausible value. Finally we adjust the zero point of the H I map, either varying the assumed zero point 
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Fig. 8. — Fits to simulated data. The y-axis gives the mean best-fit oco divided by the true value (known by construction) as a function 
of the intrinsic scatter in the dust-gas relation (x-axis). Error bars show the RMS scatter in the best- fit aco- The horizontal line at 1 
shows a perfect match between true and best-fit «co- Different colors indicate different systems. For realistic intrinsic scatter (< 0.2 dex, 
see TableO, we recover the true oco with better than 40% accuracy in all cases. 

by ±20% (in the SMC) or ±5 M© pc"^ (in the other galaxies). In detail, we vary the zero points of the IR maps 
by typically ±0.1 and 0.5 at 70 and 160//m (Itr, this depends slightly on the target). We adjust overall scale of the 
CO and H I data by ±15% (Itr) and the IR maps by ±10%. When we use a modified blackbody we pick /? randomly 
from the range 1-2 (always using a single value) and take a fraction 0.5 ±0.1 of the 70//m emission to come from 
out of equilibrium emission and so neglect it in the calculation. Each realization corresponds to a reasonable set of 
assumptions and yields a valid solution for aco a-nd (5gdr, but we emphasize that the results presented in the main 
text represent our best estimates. 

In the second test, we take Shi and Iqo ioi each data set, assume an aco: a-dd intrinsic scatter to the relation, add 
noise, and then solve for aco- We first take the real H I and CO data and assume the best-fit values of aco from Table 
[5] to be correct. We convert Jco to Eh2 and then apply the measured (Jqdr to derive a Edust for each point. At this 
point we have a plane in CO-H I-dust space that is perfectly described by a single aco and 5gdr. We add lognormal 
scatter to each quantity (equally distributed among the three axes) and finally apply noise to each axis. This yields a 
realistic approximation of real data for which we know the true aco- We do not know the intrinsic scatter a priori, so 
we try a range of values from to 0.4 dex, applying l/\/3 times this value to each axis. We repeat the exercise several 
times for each scatter, resampling the original data (allowing repeats) and re-generating the noise. 

Figure IS] shows the result of this test. The y-axis gives the best-fit aco divided by aco assumed for the simulation. 
The X-axis shows the amount of intrinsic scatter applied when making the simulated data. Each point shows the mean 
and scatter for ~ 50 tests. The statistical uncertainty reported in Table [2] is the RMS scatter (in the log) about the 
true aco for all trials with intrinsic scatter < 0.2 dex, a realistic value based on Table [2l It thus incorporates both 
the mild bias and scatter seen in Figure [H 



